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, Abstract 



The decay rates of 8 nuclides ( 85 Rr, 90 Sr, 108 Ag, 133 Ba, 137 Cs, 152 Eu, 154 Eu, and 226 Ra) were monitored by the standards group 
at the Physikalisch-Technische Bundesanstalt (PTB), Braunschweig, Germany, over the time frame June 1999 to November 2008. 
We find that the PTB measurements of the decay rate of 137 Cs show no evidence of an annual oscillation, in agreement with the 
recent report by Bellotti et al. However, power spectrum analysis of PTB measurements of a 133 Ba standard, measured in the 
same detector system, does show such evidence. This result is consistent with our finding that different nuclides have different 
sensitivities to whatever external influences are responsible for the observed periodic variations. 
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1. Introduction 

Bellotti et al. Q have recently reported the result of mea- 
surements of the activity of a 137 Cs source (Ti/2=30.08 y, 100% 
'p~ |2[]), as determined by an experiment installed deep under- 
ground in the Laboratori Nazionali del Gran Sasso (LNGS). 
They report that "no signal with amplitude larger than 9.6xl0~ 5 
at 95% C.L. has been detected," concluding that this result is "in 
clear contradiction with previous experimental results and their 
interpretation as indication of a novel field (or particle) from the 
Sun." In reviewing the case for variability, Bellotti et al. refer 
to articles by Jenkins et al. Eh, Fischbach et al. J3l, Parkhomov 
J5[0], and Javorsek II et al. Q7 ] . However, none of these articles 
cites decay rates for 137 Cs. 

Measured decay data for 137 Cs have previously been exam- 
ined for similar periodicities, and none have been found. A 
small source of 137 Cs is onboard the MESSENGER spacecraft, 
and decay data have been collected periodically from an on- 
board high purity germanium detector. These data, from just 
prior to launch on Earth to just after orbital insertion at Mer- 
cury, are consistent with no modulation of 137 Cs, as reported by 
Fischbach et al. [8J]. Ellis 09J] also reported no annual or other 
variations in measured 137 Cs decay data. What is striking about 
Ellis' result, however, is that there was an annual variation in 
the measurements of 56 Mn decay data, which were taken on the 
same detector system (over the same time period) that was used 
to measure the 137 Cs calibration standards, for which no annual 
oscillatory behavior was observed. We shall discuss the Ellis 
results in greater detail in Section [3] There is one group that 
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has reported periodicities shorter than a year in 137 Cs, | kJ UL 
but none of those experiments had a long enough duration to 
conclusively observe an annual period. 

The question of periodic or other non-random behaviors in 
nuclear decay rates has long been of interest to the scien- 
tific community Il2l - ll4ll . and the results presented in Refs. 

[I5U17I1 . as well as others, have generated renewed in- 
terest in this topic. In an effort to further explore the possible 
existence of periodicities in nuclear decays, an examination of 
historical data collected during extended studies of half-lives of 
long-lived radionuclides and of detector stability has been car- 
ried out at the Physikalisch-Technische Bundesanstalt (PTB). 
Included in this program was an analysis of 137 Cs data, as re- 
ported by Schrader llal . The goal of the present article is to 
examine the results of these PTB measurements for compari- 
son with the result of the Bellotti experiment. 

2. Analysis of PTB measurements 

Schrader iflHl reports extended half-life measurements of 8 
nuclides ( 85 Kr, 90 Sr, 108 Ag, 133 Ba, 137 Cs, 152 Eu, 154 Eu, and 
226 Ra). The 137 Cs data were collected from late 1998 to Novem- 
ber 2008. All of the measurements were made with a Ati ion- 
ization chamber (IG12/A20, Centronic 20th Century Electron- 
ics, Ltd.). In principle, these measurements could be affected 
by influences on the particular radionuclide under study, on the 
detector, or on the measuring electronics. As is clear from Fig- 
ures 2, 3, 4, and 5 of Schrader 118] , the record of measurements 
(residuals of a half-life fit) superficially resembles a scatter dia- 
gram. It follows that periodicities in the decays of any of these 
nuclides will be revealed only by some form of power-spectrum 
analysis. The purpose of the present article is to carry out such 
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an analysis for 137 Cs (residuals of a half-life fit) and also (for 
reasons that will become clear) for 133 Ba and 226 Ra. 

Since oscillations-when they occur-are typically intermit- 
tent rather than steady, it is more illuminating to examine time- 
frequency displays ("spectrograms") than simple power spectra 
I19I1 . To form spectrograms, we first prepare the data by means 
of the RONO (Rank-Order Normalization) operation [20] that 
maps the measurements onto a normal distribution, as is appro- 
priate for power-spectrum analyses such as the Lomb-Scargle 
procedure 12 lL 12211 or a likelihood procedure [23]. We then 



carry out a sequence of likelihood power-spectrum analyses of 
sections of the data. For present purposes, we have found it 
convenient to adopt sections of 500 measurements. The power, 
S , is then displayed by a color code in a time-frequency dia- 
gram. (In power-spectrum analysis, the probability of finding 
a power of S or greater at a given frequency arising from nor- 
mally distributed random noise, the null hypothesis, is given by 



[22]) 



The spectrogram formed in this way from the PTB 137 Cs data 
is shown in Figure Q] We see that there is only slight evidence 
of an annual oscillation (frequency 1 year" 1 ) between 2002 and 
2004. The feature near 0.2 year" 1 may be related to the finite 
duration of the dataset. In contrast, we show in Figure [2] the 
spectrogram formed from the 133 Ba (Ti /2 =10.551 y, 100% K- 
capture B24I0 measurements taken on the same detector system. 
This spectrogram exhibits a strong annual oscillation from 2003 
to 2005. We also see evidence of an oscillation with frequency 
close to 2 year" 1 . This could be a harmonic of the annual os- 
cillation, but it is more likely to be a Rieger oscillation (an r- 
mode oscillation with spherical harmonic indicies / = 3, m = 1), 
which is prominent in power spectra formed from Brookhaven 
National Laboratory (BNL) and PTB data I20I 
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Figure 1 : Time-frequency display (spectrogram) of measurements of the decay- 
rate of I37 Cs made at PTB over the time interval June 1999 to November 2008. 
There is only a slight suggestion of an annual oscillation from 2002 to 2003. 
The power, S , is represented by the color bar. 



We have also analyzed the PTB measurements in terms of 
"phasegrams," which are analogous to spectrograms, displaying 
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Figure 2: Spectrogram of measurements of the decay-rate of l33 Ba made at 
PTB over the time interval June 1999 to November 2008. There is evidence 
of an annual oscillation from 2003 to 2005. There is also evidence of the first 
harmonic of this oscillation. The power, S , is represented by the color bar. 



the power as a function of time and phase for an assumed annual 
oscillation. The plot derived from 133 Ba data is shown in Figure 
[3] As we expect from Figure [2] the power is found mainly 
over the time interval 2003 to 2005. The power is centered 
on a phase of approximately 0.43, corresponding to a date on or 
about 6 June. 
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Figure 3: Time-phase display (phasegram) of measurements of the decay-rate 
of 1 33 Ba made at PTB over the time interval June 1999 to November 2008. The 
phase of the annual oscillations is approximately 0.43, corresponding to a peak 
in the modulation on or about June 6. The power, S , is represented by the color 
bar. 



Comparable plots for 226 Ra (Ti /2 =1600 y, 100% a-decay 
1 25]) measurements are shown in Figures|4]and|5] We see from 
these figures that the 226 Ra results are similar to those for 133 Ba, 
but the power levels are not as strong. We note that while 226 Ra 
is 100% (T-decay, it is in equilibrium with most of its daughters, 
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several of which are p~ -decays. These fT -decaying daughters 
contribute a significant portion of the photons emanating from 
the sealed source B26I1 . Therefore, we cannot discern which 
isotope or isotopes would be the source of the observed fluctu- 
ations and note that there could be more than one. 



137 Cs standards was measured on a large, uncollimated, Nal 
system M27I1 designed for low background, high spectral reso- 
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Figure 4: Spectrogram of measurements of the decay-rate of 22s Ra made at 
PTB over the time interval June 1999 to November 2008. There is evidence of 
an annual oscillation from 2002 to 2005. The power, S , is represented by the 
color bar. 
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Figure 5: Time-phase display (phasegram) of measurements of the decay-rate 
of 226 Ra made at PTB over the time interval June 1999 to November 2008. The 
phase of the annual oscillation is approximately 0.36, corresponding to a peak 
in the modulation on or about 12 May. The power, S, is represented by the 
color bar. 



3. Discussion 

We next review in more detail previous reports concerning 
the decay rate of 137 Cs. In the Ellis experiment [9], a set of 



lution counting of human subjects undergoing in-vivo neutron 
activation analysis (NAA). The counting system, described in 
detail in Cohn et al. B27I1 . was comprised of two planar, par- 
allel 3x9 arrays of 15 cm x 5 cm matched Nal(Tl) crystals 
with low-background phototubes, for a total of 54 detectors. 
This system was located in a heavily shielded, low-background 
room at BNL. Additional measures, such as a recirculating, fil- 
tered air system, were also incorporated to reduce background, 
as described by [27]. 

The 137 Cs standard set referenced above, which was used to 
determine system stability and calibration, consisted of nine in- 
dividual sources of 0.5 fid each. Over a six month period of 
initial set-up and calibration, as reported in Cohn et al. [27], 
the standard deviation of the I37 Cs counts was 0.54%, imply- 
ing good stability, and the background over the same period 
was 2.32xl0 4 + 0.56% counts/min over a ~2.5 MeV spectrum 
window. The overall spectral resolution (ratio of peak width 
to spectrum-window width) for the 137 Cs source geometry was 
8.7%. Careful controls to monitor and quickly correct for drift 
were maintained for the system to prevent counting errors. Ad- 
ditional corrections were calculated for geometry and other pa- 
rameters which could lead to other systematic variations in the 
measured counting rate, as described in Cohn et al. [27]. The 
result was a very well designed, sensitive and stable counting 
system, with sensitivities of order 0.1 nCi. Calibration was 
performed on the system daily with the aforementioned I37 Cs 
sources. 

This counting system was utilized in conjunction with a 
broad beam neutron irradiator B28I1 . for the in vivo NAA ex- 
periments described in Ellis jjj, as mentioned above. The ir- 
radiator was specifically designed for high reproducibility j29ll . 
and comprised an array of fourteen 50 Ci encapsulated 238 Pu/Be 
sources. In 1976, a second standard set was incorporated to be 
used on a weekly basis to monitor the reproducibility of the ac- 
tivation and counting systems together [9]. This was a lucite rod 
which contained nine regions each of ~10 cm 3 volume contain- 
ing powdered manganese metal. The rod was placed in the neu- 
tron irradiator for 5 minutes, activated, and then was placed in 
the counting bed and counted for fifteen minutes. When placed 
in the counting bed, the activated regions of the rod were in the 
same locations as the 137 Cs sources with respect to the detector 
arrays 

Interestingly, the data collected on the activation product, 
56 Mn (Ti /2 =2.5789 h, 100% /T-decay, HJ), showed fluctua- 
tions with an annual period, as reported by Ellis 0], and later 
confirmed by our group in Javorsek II et al. 0]. The important 
result from Ellis' work was that even though annual fluctua- 
tions appeared in the 56 Mn data, they were not observed in the 
137 Cs standard, which was measured on the same detector sys- 
tem daily to confirm the calibration and efficiency of the detec- 
tor system. It should be noted that both the Mn and Cs standards 
were contained in similar lucite fixtures. Therefore the thermal 
expansion of the lucite, which has an expansion coefficient of 
~7.1xl0~ 5 cm/cm/°C 13111 . would have been the same for the 
two standards, and hence expansion due to temperature would 
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not be a likely source of any seasonal variations. Furthermore, 
as described in Javorsek II et al. 0], if the phase of the outdoor 
temperature is assumed to be similar to any annual variation in 
indoor temperature, it does not match the phase of the annual 
variation in the 56 Mn data reported by Ellis ]9J]. 

Based on the care and thorough effort in the design, construc- 
tion and calibration of the counting system in the Ellis exper- 
iment, aimed at reducing possible background and systematic 
variations, it seems probable that the observed differences in 
the fluctuations in the measured decay rates are intrinsic to the 
56 Mn and 137 Cs decays themselves. Of course, the possibility 
also exists that something was changed in the activation pro- 
cess, either in the production of neutrons in the PuBe sources, 
or the (n,y) capture process in the conversion of 55 Mn to 56 Mn. 
Therefore, in light of all of the examples previously mentioned 
where decay rate variations were observed in long-lived iso- 
topes, the possibility that the 56 Mn decay rate changed over the 
course of the year cannot be excluded. 

In the spirit of the above analysis it is instructive to examine 
the PTB counting system where variations in the 133 Ba count 
rate were observed, but not in 137 Cs as described in Sec. [2] As 
previously noted, the isotopes were measured on a 4n ionization 
chamber system, and the current generated in the ion chamber 
was measured with a Kiethley electrometer, model 65 17 A. This 
ionization chamber current is the result of the deposition of en- 
ergy by the interaction of the decay photons of the isotope being 
measured in the detector, generally in the detector walls. Much 
of that energy is transferred to the working gas of the detector 
(in this case argon) via Compton electrons or photo-electrons. 
Direct interaction of photons in the gas itself can also occur, but 
even at a pressure of 2 MPa (20 atm), this probability is very 
low. However, this is not true for the lower energy Compton 
photons resulting from the original wall interactions. 

One of the significant factors leading to the selection of ion- 
ization chambers in these types of measurements is their in- 
herent stability with respect to systematic and environmental 
effects. It is well known that the output of the detector is rela- 
tively insensitive to changes in the detector bias voltage, since 
the voltage is high enough to prevent recombination of the elec- 
tron/ion pairs generated, but too low to cause electron multipli- 
cation. Hence there is a fairly large width to this "plateau" such 
that the response is insensitive to voltage variations. Further- 
more, as described in detail in Jenkins et al. [32], the construc- 
tion of the detector is such that the geometry and functional pa- 
rameters of the detector will experience only negligible changes 
as a result of changes to ambient variables such as temperature, 
pressure and humidity. Possible backgrounds and their seasonal 
variations were also analyzed in Jenkins et al. rf32ll with respect 
to the annual fluctuations in the 226 Ra measurements reported 
in Siegert et al. lf33ll . These were found to be too small an effect 
to account for the observed fluctuations in the 226 Ra currents. 
It is reasonable to conclude that none of the known possible 
influences on the detector system were the likely cause of the 
fluctuations. 

It should be pointed out that the previously discussed results 
of the analysis of 226 Ra data acquired at PTB over the interval 
June 1999 to November 2008 differ from the earlier results of 



the analyses [l7ll of data acquired at PTB over the interval 
November 1983 to October 1998 Q. Notably, the decay rate 
of 226 Ra shows strong evidence of an annual oscillation in the 
earlier dataset but weaker evidence for such an oscillation in 



the later dataset (Figure 0). In fact, Schrader 11811 shows that 
strong annual periodicities were also present in 85 Kr, 108m Ag, 
152 Eu and 154 Eu measured from 1990 to 1996. 

The difference between the two sets of measurements, as de- 
scribed by Schrader 11811 . is that the current measuring system 
used in the PTB measurements prior to October 1998 was a 
Townsend balance, but for data taken after 1998, a Keithley 
electrometer was used, as described in Schrader 0~8, 34]. With 
the introduction of the electrometer, the observed annual pe- 
riodicity was reduced significantly. However, Schrader 13411 
points out in his 2007 article (8 years after the switch to the elec- 
trometer) that "for high accuracy measurements in metrology, a 
Townsend induction balance, i.e. a capacitor with voltage com- 
pensation, is used...". While it remains to be seen which current 
measurement method will prove to be the more accurate, it is 
still evident-as noted in Sec. [2]- that there clearly remains an 
annual periodicity in the 133 Ba data that is similar to those ob- 
served prior to the introduction of the electrometer, in addition 
to a weak one in the 226 Ra data, as shown in Fig. [4] 

Further support for the existence of the annually-varying pe- 
riodic behavior in the 133 Ra data is presented in Table Q] which 
lists other experiments where annual and sub-annual periodici- 
ties have been observed. What is important to note here is that 
17 of the 23 results listed in Table[T]were collected by counting 
methods utilizing pulse processing from a variety of detector 
classes, and not the current measurement from an ion chamber. 
The systematics of pulse processing are significantly different 
from the measurement of small currents, such as those of the 
PTB ionization chamber, and are not subject to the same set of 
challenges and uncertainties. 

An additional important feature of the data in Table[T]is that 
there are ten in the list that exhibit sub-annual periodicities in 
addition to the annual variations. While it is possible to attribute 
the annual periods to a "seasonal" influence with a clear annual 
variation, such as temperature, periodicities on the order of six 
months, one month, or less, are not. Furthermore, as shown in 
Refs. S [Bill El 13711 . the observed decay frequencies align 
closely with known solar periodicites. This tends to support 
the original hypotheses of Jenkins and Fischbach ifTll . Jenkins 
et al. |3], and Fischbach et al. [46], suggesting there is a solar in- 
fluence on radioactive decay rates on Earth. Even more interest- 
ing are the results presented by our group in Sturrock et al. 14311 
which indicate that the data originally presented in Parkhomov 
fl, and Parkhomov [6], exhibit a frequency structure similar to 
solar diameter measurements from the Mt. Wilson Solar Obser- 
vatory. Thus, the evidence supporting a putative solar evidence 
seems at least as reasonable as simply attributing the decay rate 
variations to unspecified environmental or systematic effects. 

Further support for concluding that the source of the vari- 
ations being observed is not simply an environmental or sys- 
tematic effect can be found in the 137 Cs and 133 Ba data pre- 
sented in Sec. [2] These measurements were made in the same 
time frame, with the same detector, and with the same current 
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Table 1 : Experiments where time-dependent decay rates have been observed. 



Isotope 


Decay Type 


Detector Type 


Radiation Measured 


Effect Observed 




Reference 


3 H 


P- 


Photodiodes 


P- 


Periodicity: 1 yr~' 




[35] 


3 H 


P- 


Liquid Scintillator 


P- 


Periodicity: 1/d, 12.1 yr -1 , 


1 yr 1 


[36] 


3 H 


P- 


Liquid Scintillator 


P- 


Periodicity: ~12.5yr _1 




[37] 


3 H 


P- 


Solid State (Si) 


P- 


Periodicity: ~2 yr -1 




[38] 




P + ,K 


Solid State (Ge) 


7 


Periodicity: 1 yr~' 




[39] 


36 C1 


P- 


Proportional 


P- 


Periodicity: 1 yr -1 , 11.7yr 


-\ 2.1 yr 1 


[3, 16,20] 


36 C1 


P- 


Geiger-Muller 


P- 


Periodicity: 1 yr -1 




[40] 


54 Mn 


K 


Scintillation 


7 


Short term decrease during 


solar flare 


[15] 


54 Mn 


K 


Scintillation 


7 


Periodicity: 1 yr~' 




[4JJ 


56 Mn 


P- 


Scintillation 


7 


Periodicity: 1 yr~' 




[9] 


60 Co 


P- 


Geiger-Muller 


P~,y 


Periodicity: 1 yr~' 






60 Co 


P- 


Scintillation 


y 


Periodicity: 1/d, 12.1 yr~' 




[42J 


85 Kr 


P- 


Ion Chamber 


7 


Periodicity: 1 yr -1 




[IB 


90 Sr/ 90 Y 


P- 


Geiger-Muller 


P- 


Periodicity: 1 yr" 1 , 11.7 yr 


-i 


[5,6,43] 




K 


Ion Chamber 


7 


Periodicity: 1 yr~' 




[18J 


133 Ba 


P- 


Ion Chamber 


7 


Periodicity: 1 yr~' 




This work 


137 Cs 


P- 


Scintillation 


7 


Periodicity: 1 d -1 , 12.1 yr" 


l 


[42] 


152 Eu 




Solid State (Ge) 




Periodicity: 1 yr~' 




[33J 


152 Eu 


P-,K 


Ion Chamber 


7 


Periodicity: 1 yr -1 




[18] 


154 Eu 


P.* 


Ion Chamber 


7 


Periodicity: 1 yr~' 




[18J 


222 Rn [ d 


a,p~ 


Scintillation 


7 


Periodicity: 1 yr" 1 , 11.7 yr 


-\2A yr 1 


[44,45J 


226 Ra W 


a,p- 


Ion Chamber 


7 


Periodicity: 1 yr -1 , 11.7 yr 


-\2.l yr 1 


[3, .17, 20] 


239p u 


P- 


Solid State 


a 


Periodicity: 1/d, 13.5 yr" , 


1 yr" 1 


[36] 



Only the count rate ratio data were available. 
Only the k photon was measured. 

Decay chain includes several primarily /^-decaying daughters which also emit photons. 



measuring electronics. It is difficult to produce a viable con- 
ventional scenario to account for this difference, particularly in 
light of the results of Ellis J9[]. Similar results were observed 
in an experiment by Alburger et al. M47I1 . who measured 32 Si 
(T 1/2 =153 y, 100% /T-decay JH) and 36 C1 (T 1/2 =301,000 y, 
98.10% ^" -decay, 1.90% K-capture fH) alternately on a dif- 
ferential gas proportional detector system f50Tl using 30 minute 
counts each for ten hours total. For each isotope, the 30-minute 
counts were aggregated into an integral count for the day. By 
taking the ratio of the two integrated counts, i.e., 32 Si/ 36 Cl, sys- 
tematic and environmental effects should have been canceled 
out. However, the ratio exhibited an annual periodicity, which 
matched, both in amplitude and phase (for the period the two 
experiments overlapped), the 226 Ra data taken at PTB IB3I1 . as 
reported in Jenkins et al. yfl. Of particular importance in the 
BNL data is that even though the two isotopes, 32 Si and 36 C1, 
were measured on the same detector on the same days, the de- 
cay data exhibited different structure, in both the amplitude and 
phase of the respective measured decay data, as discussed by 
our group in Javorsek II et al. [7] and Sturrock et al. 15111 . 



4. Conclusions 

Our purpose here has been to present new data, for 133 Ba 
and 137 Cs, which were measured on the same detector system 
for the same time period, where one isotope ( 133 Ba) exhibited a 



clear annual periodicity and the other ( 137 Cs) did not. This re- 
sult, in addition to the results of decay experiments listed in Ta- 
ble [1] indicates that the failure to observe the annual (or other) 
periodicity in one isotope does not exclude that possibility in 
others. In light of TableQ] we can state in general that our stud- 
ies to date suggest the following: (a) not all nuclides exhibit 
variability in decay constants; (b) among nuclides that do ex- 
hibit this variability, the patterns of variability (e.g., amplitude 
and phase of any oscillation) are not all the same; and (c) for 
nuclides that do exhibit variability, the patterns themselves may 
vary over time. More specifically, the results presented in Table 
Q] in support of time-varying nuclear decay rates cumulatively 
represent over sixty years of data collection, in comparison to 
the 0.5 years of data from OJJ] . 

Clearly, a new series of experiments with a variety of nu- 
clides and a variety of detectors would help to determine defini- 
tively whether the amplitudes and phases of annual oscillations 
are steady or variable and, if variable, the characteristics of the 
variability. This could lead to a determination of the cause of 
the non-random behavior exhibited by the decays of some iso- 
topes, which should eventually lead to an understanding of the 
governing mechanism. 
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